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^D , Abstract 

(N . 

We present a phase resolved spectral analysis of the X-ray binary 4U01 15+63 based upon Beppo-Sax 

f>, • satellite observations. Three strong absorption lines have been detected at all phases. We interpreted 

^-rfH ' them as a cyclotron resonant scattering. Existence of the fourth cyclotron line has been confirmed at 

^^ , some phases. The cyclotron feature is found to be strongly dependent on the phase of the source, while 

VO ' the continuum part of the spectrum depends on the phase relatively weakly. The cyclotron lines turned 

O^ I out to be very nonequidistant. The second line is almost always deeper then the first one. We discussed 

. physical conditions in the active regions of the source, under which a cyclotron feature similar to the 

(-H ' observed one can appear. 

O ■ 1 Introduction. 

C/3 , 4U0115+63 is a well-known transient massive X-ray binary source. The orbital period of the binary 

, ^/ system is 24.3 days, the eccentricity of the orbit is e = 0.34 [I]. The optical companion of the source is 

y635 Cassiopeiae, an 09 star [2]. X-ray outburst from 4U0115-I-63 has been observed first with UHURU 

,_l ! [3], and later with HEAO-1 g], CGRO/BATSE 0, RXTE [B], Ginga [7], Beppo-SAX 0. The source 

^ ' shows pulsations with the period ~3.6 s [3]. 

^^ ^ A cyclotron resonant feature in the spectrum of 4U0115+63 has been found out in 1979 9 . In 1983 

C""~~- . the second cyclotron harmonic was discovered |10) . Later the third and the forth cyclotron absorption 

^«— >> ' lines have been detected in the spectrum of the source [111 |H] . Mihara et al. fulfilled a phase-resolved 

\l^ , spectroscopy of 4U0115+63 based on Ginga observations [7]. They have found a clear dependence of 

•^^ ■ the parameters of the cyclotron feature on the phase of the source and its luminosity. 

c7^ ' By now 4U01 15+63 is the only X-ray binary having four cyclotron lines observed in its spectrum. 

OO , This unique property make the source 4U0115+63 a touchstone for theories of the gyromagnetic emission 

^^ ' of X-ray pulsars. If only one weak line is visible in the spectrum of a pulsar, usually it can be fitted well 

^ I enough by a very wide group of theoretical models. A big number of feature's parameters in the case of 

4U0115+63 could give as an opportunity to choose from this abundance of models the appropriate one. 
Here we report on the phase-resolved spectral analysis of 4U01 15+63 based on Beppo-Sax satellite 



X 



JH , observations of the source's outburst in 1999 March. We concentrated particulary on the cyclotron 

feature and its evolution with the phase. Obtained results allowed us to analyze the conditions in the 
active region of the source and physical processes leading to the formation of the cyclotron lines. 

2 Observations and spectral analysis. 

2.1 Time analysis. 

This work is based upon the observations performed by Beppo-Sax satellite on 1999 March, 26 (OP6714) . 
Observation time was equal 17''12.5''. Since the first cyclotron harmonic in the spectrum of 4U01 15+63 
is near ~12 keV [101 j^ we used the data from the spectral interval from 2 keV to 75 keV. This energy 
range is covered by three Beppo-Sax instruments: MECS2(2-11 keV), HPGSPC (7-30 keV), PDS(15- 
75 keV). More detail information about these instruments can be found in [121 1131 [T4l [8]. 
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Figure 1: Phase profiles of tlie MECS2 data for the energy ranges 4 — 11 keV. 
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Figure 2: Phase profiles of the HPGSPC data for the energy ranges 7 — 30 keV. 
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Figure 3: Phase profiles of the PDS data for the energy ranges 15 — 75 keV. 
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Figure 4: Phase profiles of the source 4U01 15+63 for the energy ranges 2 — 6 kcV. 
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Figure 5: Pulse profile for the energy range 4—11 keV. Phase region selection is shown. 



In order to obtain the energy-resolved pulse-profiles and phase-resolved energy spectra the photon 
arrival times were first corrected to the Solar System barycenter. By using period folding [15J, we have 
found the period of the pulsar to be equal 3.6143 ± 0.0002 s. No variations of the period of the source 
have been detected. Obtained pulse profiles for different energy ranges are represented on Fie: IH3l A 
main soft peak (in the phase interval 0,2-0,75) with maximum at 0,4 and a hard peak (in the phase 
interval 0,75-1,2) with maximum at 1,05 are plainly distinguishable. Hard peak is relatively narrow, 
and its shape remains almost constant over a wide range of energies (above 4 keV). On the contrary, 
the form of the main soft peak depends highly from energy; over 30 keV it almost disappears. 

To examine the situation in more details we plotted the softest part of the emission with higher 
energy resolution (Fig|4)| .On this figure a new feature appears: a secondary soft peak with maximum 
at the phase 0,85, while the hard peak disappears. Moreover, the shape of the main soft peak changes 
drastically. On the energies above 3 keV it has a clear maximum at the phase 0,4. and descending 
region between 0,4 and 0,7. But over the energies 2-3 keV it becomes quite flat without any clear 
maximum. It might be accounted for the influence of the above-mentioned secondary soft peak, but if 
we consider phase profile evolution from 5 to 2 keV, it turns out, that intensity difference between the 
maximum value at the phase 0,85 and the minimum value at ~0,75 remains almost the same, while 
intensity difference between the minimum value at 0,75 and the value at 0,7 rises steeply. It means 
that the appearance of the secondary soft peak and the change of the shape of the main soft peak at 
low energies are two different phenomena. 
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Table 1: Continuum properties in the phase interval 0,4-0,8 
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Figure 6: The spectrum of the pulsar at the phase 0,55, in the soft peak. 



2.2 Spectral analysis. 

In order to perform a phase-resolved analysis of the spectrum of the pulsar, all phase interval has been 
divided on ten equal phase regions. The phase region selection is represented on FigO 

The source 4U0115-I-63 has a power- law spectrum with high-energy cutoff; this continuum is strongly 
distorted by cyclotron absorption |10|[5]. It turned out necessary during our analysis to divide the whole 
phase interval on two regions and to use different continuum models for them. First region covers phase 
interval from 0,4 up to 0,8 and approximately corresponds to the descending edge of a lower peak of 
the pulse-profile; second region covers all other phases. For the first region it was enough to use a simple 
power-law with exponential cutoff model 



/(£)(xi5%xp(-^) 



(1) 



For the second region it turned out to be insufficient. The continuum spectra within the main peak 
have been fitted by the following model: 



(2) 



f(E)^^ B£«2exp(-|), E>E^^, 

Numerical coefficients A and B have been chosen so that to make the function continuous. A fracture 
on the energy Ecut has been smoothed by a polynomial function. 

In order to represent the cyclotron feature these continuum models have been multiplied on the 



so-called Gaussian filters, one for each line. 



G„(£) = l-D„exp 
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(3) 



data and folded mode 
1ECS2_Energy_f.fits HPGSPC_Energy_f.fits PDS_Energy_f.fits 
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Figure 7: The spectrum of the pulsar at the phase 1,05, in the center of the hard peak. 

Here n is the number of harmonic, Ei is an energy of the first harmonic. Parameters A'^ define positions of 
the other harmonics. Initially they were taken equal to n, but they were considered as free parameters in 
order to represent unequidistant cyclotron lines. Three cyclotron lines have been detected at all phases. 
The fourth harmonic turned out to be necessary only at the phases 0,9-1,2. No iron line at 5-7 keV 
has been detected. 

Described model turned out to be as a whole quite acceptable to fit the spectra, though in some phase 
intervals certain problems appeared (see Discussion). Results, obtained with the model, are represented 
on tables [l][5] and figures 1 81131 



3 Results and discussion. 

3.1 Energy-resolved analysis of the phase profiles. 

As it has been described in the previous section, pulse profile of the source contains two main features: 
the hard peak and the main soft peak. This brings up a question: Are these features produced by 
different magnetic poles of the neutron star, or they are a pencil and a fan beam from the same pole. 
Let us discuss both these hypotheses. But at first we remark, that if the emission has a fan directional 
diagram, an observer should cross its plane twice a period. Therefore we can expect a fan beam should 
produce two features on a pulse profile. And so the hard peak is the most probably a pencil beam (there 
is only one hard peak in the pulse profile), while the two peaks could either be produced by the same 
fan beam (the phase distance between them is almost equal 0,5), or be two independent components, 
produced by different magnetic poles. 

If the hard peak and the main soft peak are formed in the same active region (this supposition 
is called hereafter hypothesis 1), then the main soft peak is, of cause, a fan beam (the phase interval 
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Figure 8: Power-law indexes. Detail description of the parameters a,ai and 02 see on the formulas ([T]) and 
0. 



between it and the hard peak is big enough), and the secondary soft peak is almost indubitably the 
second impulse of this fan. The hard and the main soft peaks together cover nearly the whole period of 
the pulsar. It means that one active zone of the neutron star is visible all the period, while an emission 
of the other does not reach an observer at all. This situation is possible only if the angle between our 
line of sight and the axis of rotation of the neutron star is small enough. 

The hard peak is situated strongly asymmetrically with respect to the two soft peaks, which is to 
say that the pencil beam is inclined noticeably with respect to the axis of the fan beam. Moreover, the 
spectral properties of the main and the secondary soft peaks are quite different, i.e spectral characteristics 
of the fan beam depend highly from azimuth. Surely, it is possible only in the accretion occurs non- 
axis-symmetrically. It means, that the magnetic field of the pulsar is greatly non-dipole. 

Let us now consider the situation, when the hard peak and the main soft peak are formed in different 
active regions (hereafter hypothesis 2). Then the main soft peak and the secondary soft peak could be 
components of a fan beam forming on the second magnetic pole of the neutron star {hypothesis 2a). 
In conclusion, the situation is possible when the main soft peak is formed by a pencil beam from the 
second magnetic pole, while the secondary soft peak is only a softer and wider component of the hard 
peak {hypothesis 2b). But in any case there are strong arguments that the field of the pulsar differs 
drastically from the dipole one. The hard peak is situated strongly asymmetrically with respect to the 
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Figure 9: Temperature Ti and T2, and cutoff energy -Ecwf 
Ecut see on the formulas (H]) and ([2])- 



Detail description of the parameters Ti ,r2 and 



two soft peaks (this argument is valid if we assume hypothesis 2a) and the phase interval between the 
hard peak and the main soft peak is well below 0.5 (if we assume hypothesis 2b). 

Under consideration of the hypothesis 2 a question appears: why we can not see the hard component 
from the second magnetic pole? There are two probable explanations of this fact. As the hard component 
is more concentrated, it is possible that the geometry of the source is so that a hard component from 
the second hot spot on the neutron star's surface (corresponding to the soft peak) exists but is invisible 
for us. But it is also possible, that emission of this hot spot doesn't contain a hard component at all. As 
it has been already mentioned, magnetic field of the pulsar differs strongly from the dipole one. Under 
such conditions the physical parameters of the accreting streams and, consequently, of the emitting 
zones are absolutely different for the magnetic poles of the neutron star. It can led to the absence of a 
hard component from one of the poles. 

Both hypotheses {hypothesis 1 and hypothesis 2) have arguments for and against. If we assume 
that the hard peak and the main soft peak are formed on different active regions, then (considering the 
above-mentioned departure of the magnetic field from the dipole one) we might expect that the cyclotron 
features in the spectra of the hard and the main soft peaks have difi'erent fundamental energies. As seen 
from Fig lfOI this is not the case. Changes of the first cyclotron line's energy doesn't correlate with the 
peaks of emission. This testifies in favour of formation of all the emission in the same region {hypothesis 
1). But the hypothesis 1 faces with an another trouble. As it will be discussed below, the emission 
of the pulsar likely appears as a result of comptonization of the emission by hot plasma. Spectral 
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Ta 


0,85 


-0,164±0,16 


-2,65±0,38 


12,0±0,38 


25,0±5,5 


0,95 


0,0525±0,027 


-0,349±0,043 


6,78±0,10 


10,8±0,24 


1,05 


0,539±0,20 


-0,966±0,85 


9,00±0,19 


12,9±4,1 


1,15 


0,256±0,023 


-1,33±0,052 


8,69±0,10 


17,7±0,61 


1,25 


-0,0926±0,033 


-1,53±0,10 


8,40±0,21 


18,2±1,1 


1,35 


-0,219±0,065 


-1,68±0,18 


7,91±0,33 


16,0±1,6 



Table 2: Continuum properties in the phase interval 0,8-1,4 

properties of the obtained emission depends upon the temperature of the plasma [16]. The secondary 
peak disappears already at 3-4 keV, the main soft peak - at 25-30 keV; The hard peak is visible up 
to the highest energies. If we assume the hypothesis 1, it means that in a single active zone there are 
tree regions of plasma with absolutely different temperatures, which on the one hand have an optical 
depths great enough to produce remarkable features on the pulse profile (and each feature has their own 
directional characteristic), but on the other hand each region does not infiuence at all on the emission 
of the others. This situation is, of cause, possible, but it seems a bit artificial. It would be more natural 
to suppose, that the emission is formed in two different zones with different physical conditions. 

3.2 Spectral analysis. 

The source 4U0115+63 has a power-law spectrum with high-energy cutoff . This shape of spectrum as 
a whole is typical for the X-ray binaries and likely appears as a result of comptonization of the emission 
by the hot plasma of the accretion column. We tried to fit the continuum by the spectra suggested in 
[16] . These spectra represent results of comptonization by a simple-geometry plasma slab. Our attempts 
failed. It means that the geometry and the temperature distribution in the active region of the source 
are rather complicate and they probably consist of several layers of plasma with various temperatures 
and optical depths. 

The cyclotron lines exist definitely at all the phases of the pulsar, but parameters of the spectrum 
are very phase-dependable. The parameters are represented on tables [l][2] First, second and third 
harmonics have been detected at all the phases of the source; fourth harmonic was visible only in the 
hard peak. 

The cyclotron feature is the most pronounced not on the maxima of luminosity but on the descending 
edges. Moreover, the fundamental cyclotron energy also has a maximum on these edges. One of the 
possible ways of explanation of this effect is the following: an emission propagating exactly along the 
accretion column interacts with plasma on a very large interval of heights. The main energy of a 
cyclotron line is formed well above the surface of the neutron star in a region with respectively weak 
magnetic field. On the contrary an emission propagating aslant with reference to the column leaves the 
column respectively deep, but its path in these deeper layers is longer. Consequently, we can expect 
that in this case that cyclotron line would be deeper but narrower and its main energy higher. One 
can see that it is really so: the widths of the lines have local maxima correlating with the maximum of 
luminosity instead of the descending edge. But then a question appears: why this effect acts only on the 
descending edge and not valid on the ascending one? The cause of this could be the above-mentioned 
asymmetry of the accretion column appearing as a result of out-of-dipole magnetic field of the pulsar. 
The difference between the spectra of descending and ascending edges can be considered in its turn as 
an addition evidence of a departure of the field from a dipole. 

The fundamental cyclotron energy is not stable: its variations are of the order of 10% and they are 
near correlating with the variations of luminosity. Harmonics are very unequidistant. This effect is the 
most distinct on the second harmonic. The ratio between energies of the second and the first harmonics 
is as small as 1.6 at some phases. 

The situation with the theoretical interpretation of these effects is not so simple. The effects are 
usually interpreted as a result of complex geometry of the emission forming region and magnetic field 
variations in it. Let as analyze this possibility. 

The neutron star is considered to have a radius ~10 km and a dipole magnetic field (for the estimating 
consideration this supposition is quite acceptable). The luminosity of the source 4U0115-I-63 is ^ 
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0,43±0,031 
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0,53±0,055 


0,58±0,063 


— 


0,65 


0,28±0,25 


0,32±0,10 


0,51±0,070 


— 


0,75 


0,28±0,095 
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— 


0,85 
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0,46±0,061 
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0,048±0,0090 


0,44±0,017 


0,27±0,041 


0,32±0,025 


1,05 


0,47±0,14 


0,59±0,087 


0,35±0,043 


0,43±0,25 
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0,23±0,0095 


0,50±0,0071 


0,48±0,37 


0,26±0,033 


1,25 


0,21±0,020 


0,46±0,010 


1.0±3,6 


— 


1,35 


0,23±0,026 


0,37±0,016 


0,55±0,37 


— 



Table 3: Depths of the cyclotron lines. 
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0,45 


2,77±0,21 


4,41±2,1 


6,41±5,8 
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0,55 


3,96±1,1 


4,41±3,0 


6,41±6,9 
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3,50±1,4 


4,41±8,6 


6,41±7,4 
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3,16±0,58 


4,41±3,4 


6,41±4,8 
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3,72±0,66 


3,97±0,51 


4,23±1,4 
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0,749±0,22 


3,21±0,14 


1,86±0,60 


8,98±1,9 
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4,41±4,04 


3,69±0,68 


2,12±0,58 


12,7±5,1 


1,15 


1,77±0,13 


2,02±0,065 


1,09±0,97 


7,69±1,9 


1,25 


1,88±0,20 


2,19±0,099 


0,429±0,77 


— 


1,35 


1,62±0,11 


2,68±0,19 


0,805±0,64 


— 



Table 4: Width of the cyclotron lines. 
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0,45 


12,2±0,10 


1,72±0,038 


2,55±0,092 
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0,55 


12,56±0,28 


1,69±0,052 


2,55±0,11 
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0,65 


12,39±1,2 


1,51±0,056 
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0,75 


12,19±0,44 


1,62±0,034 


2,84±0,12 
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11,15±0,17 


1,77±0,020 


3,08±0,089 
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0,95 
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1,9±0,027 
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3,77±0,11 


1,05 


12,23±0,76 


1,84±0,12 


2,69±0,17 
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12,8±0,73 


1,86±0,0099 


2,83±0,035 


4,16±0,10 


1,25 


12,11±0,13 


1,99±0,020 


3,04±0,053 


— 


1,35 
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2,04±0,015 


2,77±0,034 
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Table 5: Fundamental cyclotron energy and positions of the cyclotron lines. 
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Figure 10: Energy of the first cyclotron harmonic Ei 



W"^^ erg/sec. When determined from the cyclotron energy by the usual nonrelativistic formula 



Ei = n 



eH 

TUeC 
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(4) 



the magnetic field of the source appears to be 1.1 ■ 10 Gs. There is a well-known formula for the 
angular size x of the hot spot on the X-ray pulsar's surface if its magnetic field is a dipole. 



X~sinx = 7.5- 10"''// 7L7 



(5) 



In the considering case we obtain x — 0.06. Relative variation of the magnetic field (and hence of the 
cyclotron energy) within the boundaries of such small a spot is of the order of -^ ~ 1.4 • 10~^. Of 
course it is enough to explain neither widths of lines nor their arranging. Consequently all these effects 
should be considered as a result of variations of the magnetic field along the height of the accretion 
column. 

In the average the second harmonic is the strongest. The first harmonic is the most pronounced 
only in the region out of the major peak (phases 0,45-0,85). At the phase 1,15 where the cyclotron 
feature is the most distinct the second and the third harmonics have almost the same depths which 
is significantly higher then the depths of the first and the fourth one. These ratios between depths 
of different harmonics looks strangely. As it was shown in TS", cyclotron cross sections averaged over 
directions and polarizations under considering physical conditions are related as 
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Figure 11: Depths of the cyclotron hnes. 

Even for the direction perpendicular to the magnetic field this ratio is 

(Ti : (72 : erg = 200 : 18 : 3 



(7) 



One can see that in any case the first harmonic should be significantly stronger then the others. Mean- 
while at the main peak the second and the third harmonics are stronger then the first one. The most 
believable explanation of this fact was given in the series of articles (see [T5], [20]) [S]: [12], [23] )• The 
authors explain the contradiction by the influence of the Raman scattering. Under the physical condi- 
tions existing in the accretion column the probability of the transition of an electron from the second 
Landau level through the branch (2 ^ 1 ^ 0) is approximately seven times higher [24) then through 
the branch (2^0). An influence of the true absorption is negligible [25]. As a result a photon of the 
second harmonic being absorbed decays on two photons of the flrst one. This spawning leads to the 
signiflcant decreasing of the first harmonic's depth. 

But this mechanism could be acceptable only for the second harmonic: for the ratio between third 
and second harmonics it fails. In fact, the integrated probability of the branches (3 ^ 1 ^ 0) and 
(3 — > 2 — > 0) is only 29%, and only one photon of the second harmonic appears instead of two photons 
in the previously described case. Furthermore even this rare photon has no chance to reach the distant 
observer. As it was shown in [21] . Raman scattering of photons of the third harmonic can infiuence 
effectively on the spectra of previous harmonics only if an optical depth of the scattering electrons on 
the third harmonic is at least comparable with 1. But in this case their optical depth on the second 
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Figure 12: Widths of the cyclotron hnes. 



harmonic should be ~6. The spawned photon of the second harmonic is sure to experience the second 
scattering and decay on two photons of the first harmonic. So the efficiency of the spawning of photons 
of harmonics of the orders higher then first is negligible. But ratios of their depths in the observing 
spectrum are also very far from the theoretical relations between cyclotron cross-sections. 

The nonequidistancy of the harmonics is also very difficult for the theoretical explanation. Owing to 
the above-mentioned relations between the cyclotron cross-sections the magnetized substance is much 
more transparent on the energy of the second harmonic then on the energy of the first one. Consequently 
the second harmonic should appear in the deeper regions of the accretion column with the stronger 
magnetic field. We might expect the distance between first and second harmonics to be bigger then the 
harmonic one. Meantime the distance is smaller than the harmonic: the ratio between energies of the 
second and the first harmonics sometimes is as small as 1.6. The same result have been already reported 

■mm 

The third harmonic is grotesquely deep at the phase 1,35. Moreover the confidence intervals of 
the D3 at phases 1,25-1,45 are too big. The reasons of this effect is that the real line shapes differ 
distinctly from the Gaussian. The same result has been already reported 0. Unfortunately there is 
now no universally accepted view of what shape should have a cyclotron line in a spectrum of an X-ray 
pulsar. A Gaussian model is one of the most commonly-using, mainly because of its simplicity. We 
tried different kinds of models for the lines and the Gaussian turned out to be the best. Moreover, the 
emission at the phase 1,25 is a mixture of emissions of the hard peak and the main soft peak and it 
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Figure 13: Positions of the cyclotron lines. See the significance of the parameters N2 and N^ in the formula 



m 



makes its spectrum very complex. One of the possible ways to avoid this difficulties is to suppose that 
the cyclotron feature is formed by cyclotron radiation of ultrarelativistic, anisotropic electrons |17) . 

4 Conclusion. 

To summarize; We studied the 4-75 keV pulse phase resolved spectrum of 4U0115-I-63 observed by 
Beppo-Sax. Analysis of it suggests that the magnetic field of the neutron star differs drastically from 
the dipole one. We confirm the strong dependence of the parameters of the cyclotron feature in the 
spectrum on the phase, while the continuum part changes rather weakly. We also confirm the existence 
of the fourth cyclotron line near ~ 45 keV. Gaussian shape is not quite good to fit the lines; further 
investigations are needed to propose some better model. 

Consideration of the physical processes of the cyclotron feature formation points to some serious 
difficulties of the absorption model. In particular, it is unclear how to explain the observing ratios 
between depths and central energies of cyclotron lines in the context of it. Thus the generally used 
absorption model of the cyclotron feature in the spectrum of the source can not be considered as an 
indisputable. 
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